Variational method for lattice spectroscopy with ghosts 
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We discuss the variational method used in lattice spectroscopy calculations. In particular we 
address the role of ghost contributions which appear in quenched or partially quenched simulations 
and have a non-standard euclidean time dependence. We show that the ghosts can be separated 
from the physical states. Our result is illustrated with numerical data for the scalar meson. 
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QUENCHED CORRELATORS AND THEIR 
SPECTRAL REPRESENTATION 



Spectroscopy of excited states is still a quite challeng- 
ing problem for lattice QCD. The reason is that euclidean 
two-point functions have a sum of infinitely many ex- 
ponentials in their spectral representation and it is a 
non-trivial task to extract the subleading terms corre- 
sponding to the excited states. For quenched or partially 
quenched calculations the situation is complicated fur- 
ther by additional unphysical contributions: When ig- 
noring the fermion determinant, the rj' also becomes a 
Goldstone boson, and this massless flavor singlet state 
can couple in various processes through hairpin diagrams. 
The corresponding effects were studied in (jy| for the 
quenched scalar propagator. It was demonstrated that 
the T]' leads to additional contributions to the scalar prop- 
agator which have a negative sign and a non-standard 
f-dependence. Such terms, often referred to as "ghosts" , 
compete with the exponential decay coming from the 
scalar, and make the extraction of the scalar mass cum- 
bersome. For the nucleon system the effects of ghost 
states and the problems they cause in the extraction of 
the Roper resonance were discussed in . 

In principle, in a quenched calculation, correlators for 
all quantum numbers may be infested by ghosts, possibly 
coming from bound or scattering states with 77', or states 
involving even several rj'. Thus, one has to find tools 
for reliably removing these contributions, in particular, 
if one is interested in excited states, where the influence 
of the ghost states is even more devastating. 

A powerful tool, developed for the analysis of excited 
states, is the variational method The main idea is 

to work with several linearly independent interpolators 
Oi, i = 1, 2, . . . r, all with the quantum numbers of the 
desired state, and to compute the cross-correlations 



In full euclidean lattice QCD the correlators have the 
spectral decomposition 
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where the coefficients are given by 



= (O|0.|n) 
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In Eq. 121) Vn denotes the energy of the state \n). 

In the presence of ghosts we must augment the decom- 
position by additional terms which take into account 
the possible negative sign and the different i-dependence 
of the ghost contributions. Since, in principle, there are 
infinitely many ghost states, their contribution also gives 
rise to an infinite tower of states 
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Here Wn is the energy of the n-th ghost state and the 
function f'^'^\t) describes the non-exponential part of its 
euclidean time dependence, taking into account also a 
possible negative sign. We stress that for our analysis we 
do not have to know the precise form of /^"^ (t). Certainly 
it is not simple to present a form which encompasses all 
possibilities that might occur. However, to give an exam- 
ple, we remark that the linear behavior f(t) = —{a + ht) 
was suggested for the system considered in 2]. 

In order to combine the contribution (0J from the 
ghosts and the original decomposition we represent 
our correlation matrix as 
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where 



C{t), 



(O,(i)O,(0)) , i,j =1,2, 



(1) 



_ff(")(t) — 1, En — Vnp for proper states , 
i7(")(t) = /("^'(t) , = W„G for ghosts. (6) 
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Here we have ordered the energy levels such that (assum- 
ing non-degenerate values for the E'^"'>) 



< El < E2 < E3 .. 



(7) 



implying that the ghost state energies are interleaved 
with the energies of proper states. A typical situation 
would be 



< Vi < Wi < V2 ... , 



(8) 



where the energy Wi = E2 of the first ghost state comes 
after the ground state energy Vi = Ei, but is below the 
first excited energy level V2 = E3. In such a case the 
ghost contribution overlays the exponential decay from 
V2 and makes the extraction of this number from a single 
correlator quite difficult. 



II. VARIATIONAL METHOD WITH GHOST 
CONTRIBUTIONS 



are given by JH)). For linearly independent interpolators 
Oi, that couple to the first r states, full rank follows 
from the fact that the states |n) are orthonormal. For 
the quenched approximation one expects that the ghost 
states are linearly independent of the proper states such 
that again for linearly independent Oi the matrix a["^ 
has full rank. 

Inserting (|12|l in one finds that for all i 
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with the coefficients p^"'*^) given by 
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From the full rank of a] it follows that the vectors 
r, are linearly independent, implying 
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Having generalized the formula for the spectral repre- 
sentation from the original form |2l to the ansatz 
which is capable of describing also the ghost contribu- 
tions, we can now proceed as in 0. In particular, we 
consider the generalized eigenvalue problem 
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On the right-hand side we use the correlation matrix at 
some fixed < t Xo normalize the eigenvalue problem. 
We now show that the eigenvalues \{t) of the generalized 
eigenvalue problem are given by 



-(t-to)Vfc 
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l + C'(e-(*-*'')'^'') (10) 



for proper states. 



for ghosts. 



where E^ is the energy of the fc-th state (proper state or 
ghost) and = E^+i — E). is the difference to the mass 
of the r + l-st state, r being the number of interpolators 
used for the correlators. 

The proof proceeds by considering the generalized 
eigenvalue problem 
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(15) 

Using the linear independence of the eigenvectors tp^'^^ 

(n) 

and the full rank of Uj , it follows that for each n there 
exists a k, such that p^"'*"') is non- vanishing. Conse- 
quently, Eq. lfT31) implies that the eigenvalues of the re- 
duced problem Hll() are given exactly by the leading term 
in (^ni)- The full matrix C{t) is obtained from C{t) by 
adding terms of 0{exp{—t Er+i)) which can be treated 
perturbatively 0, giving rise to the correction in (|1(J|I . 

Before we come to discussing an application of the gen- 
eralized variational method, let us briefly address once 
more the implications of the result (|10|l : In the leading 
term each eigenvalue couples to only one state (ghost 
or proper state). Thus, ghosts can be cleanly disentan- 
gled from the proper states, up to the correction term 
C'(exp(— t Afc)). We stress that the form of the euclidean 
time dependence for the ghosts, /*•''■'(<), needs not be 
known or modelled for the application of the method. In 
the example we discuss below, the ghost state is easily 
identified by its non-standard ^-dependence. 

We remark that additional information for the identifi- 
cation of the ghost states comes from the standard eigen- 
value problem, i.e., when the normalization with C(to) is 
omitted on the right-hand side of JHJl . Following the argu- 
ments given in the appendix of 4] , one finds that for the 
standard eigenvalue problem the eigenvalues are given by 



for the hermitian r x r matrix C(<), obtained by truncat- 
ing the spectral sum after the r-th term. 
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We assume that the coefficients a|"\ i,n = 1,2 
form a matrix of full rank. For the proper states the a. 
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r, 
(«) 
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Since here the factor H^''\t) is not normalized by 
H'^'^Hto), the sign information is not lost and an over- 
all sign indicates a ghost. However, this information is 
obtained at the cost of a much larger correction term, 
since in (|16ll Sk is the distance to the nearest energy level 
which is usually much smaller than appearing in H1U|I . 
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FIG. 1: Euclidean time dependence for the diagonal entries of the correlation matrix (top row) and for the eigenvalues of the 
generalized eigenvalue problem (bottom row). The quark masses are am — 0.02, 0.04 and 0.06 (left to right). 



III. ILLUSTRATION OF THE RESULT FOR 
THE SCALAR MESON PROPAGATOR 

We illustrate our result using the scalar meson as an 
example. In particular, we work with the three interpo- 
lators 



"^n dn , tin du 



(17) 



The subscript n denotes a narrow quark source (or sink), 
while w is used for wide. The narrow and wide quark 
sources (sinks) are constructed via different amount of 
Jacobi smearing |^] . Such interpolators with mixed width 
for the quark sources have been utilized in lattice spec- 
troscopy of excited mesons and baryons and details of the 
source preparation can be found in 

We remark, that the interpolators used here are not 
meant to resolve the issues currently discussed for the 
scalar meson (see e.g. 0,0). They only serve to illustrate 
our result for the system where the importance of ghost 
contributions was first discussed Q. A detailed analysis 



of the spectroscopy results for the interpolators H17() will 
be presented elsewhere. 

Using the three interpolators (|17|l we work with a 3 x 3 
correlation matrix C{t). We evaluate it on 100 quenched 
gauge configurations on a 20^ x 32 lattice using the chi- 
rally improved lattice Dirac operator j^. The gauge ac- 
tion is the Liischer-Weisz action at f3 = 8.15, correspond- 
ing to a lattice spacing of a = 0.119 fm determined from 
the Sommer parameter We present results at three 
different quark masses am = 0.02, 0.04, 0.06, gi ving rise 
to pion masses of ~ 400, 550 and 670 MeV [ll||. All 
errors given are statistical errors determined with single 
elimination jackknife. 

In the top row of plots in Fig. ^ we show the euclidean 
time dependence for the diagonal entries of the correla- 
tion matrix. We remark that the off-diagonal elements 
show qualitatively the same behavior, but we omit them 
in the plot to avoid overcrowding it. For the smallest 
quark mass the entries of the correlation matrix become 
negative at t/a = 3 or 4, showing the strong influence 
of the ghost contribution. For am = 0.04 the effect is 
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still visible, but the correlators become negative only at 
values t/a — & or 7. At the largest quark mass the cor- 
relators do not become significantly negative, i.e., within 
error bars are compatible with zero for t/a > 10. This 
indicates that the role of the ghost state is much weaker 
at larger quark masses. This is as expected from the 
fact, that the mass of the would-be Goldstone particles, 
and thus of the ghost states, increases very fast with the 
quark mass. 

The bottom row plots show the time dependence for 
the three eigenvalues of the generalized eigenvalue prob- 
lem. They behave as expected from formula (|10|l . Only 
one of them, \'^^\ shows the behavior characteristic of 
ghosts. The other two eigenvalues A*-^-* and A^^^ remain 
positive and are compatible with a single exponential 
decay. The sign change of A^^^ is most prominent at 
am = 0.02 and the overall time dependence is similar to 
the behavior of the corresponding diagonal elements of 
the correlation matrix in the top left plot. For am — 0.04 
the ghost-like behavior of A^'^-' is still visible, while for 



am = 0.06 it is no longer significant. 

Comparing the top and bottom row plots of Fig. 1 il- 
lustrates the strength of the method: For the two masses 
am = 0.02 and 0.04, where the single correlators show a 
clear ghost contribution, the ghost appears in only one of 
the eigenvalues, while the other two are consistent with a 
dominant single exponential behavior. Our result shows 
that with the variational method ghost and proper states 
each dominate individual eigenvalues and the masses of 
the physical states can be fit without having to know or 
model the form of the ghost contribution. 
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